Errors in trapped-ion quantum gates due to spontaneous photon scattering 
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We analyze the error in trapped-ion, hyperfine qubit, quantum gates due to spontaneous scattering 
of photons from the gate laser beams. We investigate single-qubit rotations that are based on stim- 
ulated Raman transitions and two-qubit entangling phase-gates that are based on spin-dependent 
optical dipole forces. This error is compared between different ion species currently being investi- 
gated as possible quantum information carriers. For both gate types we show that with attainable 
laser powers the scattering error can be reduced to below current estimates of the fault-tolerance 
error threshold. 
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I. INTRODUCTION 

Quantum bits, or qubits, that are encoded into internal 
states of trapped ions are an interesting system for Quan- 
tum Information Processing (QIP) studies 0,0]. Internal 
states of trapped ions can be well isolated from the en- 
vironment, and very long coherence times are possible 
0, H, 0) Hi • The internal states of several ion-qubits can 
be deterministically entangled, and quantum gate s can be 
carried out between two ion- qubits @, 0, n MM, nn m • 

Among different choices of internal states, qubits that 
are encoded into pairs of ground state hyperfine or Zee- 
man states benefit from negligible spontaneous decay 
rates [2J. The small energy separation between the two 
states of such qubits (typically in the radio-frequency or 
microwave domain) allows for phase coherence between a 
local oscillator and a qubit superposition state over rela- 
tively long times 0, S El Q3] ■ 

The quantum gates that are performed on hyperfine 
ion-qubits typically use laser beams. Since light couples 
only very weakly to the electron spin, spin manipulations 
rely instead on the spin-orbit coupling of levels that are 
typically excited non resonantly through allowed electric 
dipole transitions. Spin manipulations therefore require 
a finite amplitude in the excited electronic state, and 
spontaneous scattering of photons from the laser beams 
during the gate is inevitable. 

Fault-tolerant quantum computation demands that 
the error in a single gate is below a certain threshold. 
Current estimates of the fault-tolerance error threshold 
range between 10 -2 to 10~ 4 [l5|, [TH, [13] • These estimates 
rely on specific noise models and error-correction proto- 
cols and should be considered as guidelines only. How- 
ever, the general view is that for fault tolerance to be 
practical, the error probability in quantum gates should 
be at least as small as 10 -4 . It is, therefore, worth ex- 
ploring the limitations to the fidelity of quantum gates 
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performed on trapped ions with laser light, using this 
level of error as a guideline [18| ■ 

In his 1975 paper [r|, Mollow showed that the effect of 
a quantum coherent field on an atom is equivalent to that 
of a classical field plus a quantum vacuum field. The error 
due to the interaction with light can be categorized into 
two parts. The first is the error due to noise in classical 
laser parameters, such as intensity or phase [2t \20\. The 
second part originates from the quantum nature of the 
electromagnetic field and is due to vacuum fluctuations, 
i.e. the spontaneous scattering of photons pll. [2^|. 



A. Ion-qubit levels and transitions 

Most ion species considered for QIP studies have a sin- 
gle valence electron, with a 2 Si/ 2 electronic ground state, 
and 2 Pi/2 and 2 P3/2 electronic excited states. Some of 
the ions also have D levels with lower energy than those 
of the excited state P levels. Ions with a non zero nu- 
clear spin also have hyperfine structure in all of these lev- 
els. A small magnetic field is typically applied to remove 
the degeneracy between different Zeeman levels. Here we 
consider qubits that are encoded into a pair of hyperfine 
levels of the 2 Si/2 manifold. Figure [T] illustrates a typical 
energy level structure. 

To allow for a straightforward comparison between dif- 
ferent ion species, we investigate qubits that are based 
on clock transitions, i.e. a transition between the \F — 
I - 1/2, m F = 0) = | T) and |F = 1+1/2, m F = 0) = | 1) 
hyperfine levels in the S1/2 manifold of ions with a half 
odd-integer nuclear spin /, at a small magnetic field. For 
this transition, the total Raman and Rayleigh photon 
scattering rates, as well as the Rabi frequency, are inde- 
pendent of /, and the comparison between different ion 
species depends only on other atomic constants. Super- 
positions encoded into these states are also more resilient 
against magnetic field noise 0, [f| [ll[ . Even though our 
quantitative results apply only to this configuration, for 
other choices of hyperfine or Zeeman qubit states (in- 
cluding those with 1 = 0) the results will not change 
significantly. 
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FIG. 1: Relevant energy levels (not to scale) in an ion-qubit, 
with nuclear spin I. The P1/2 and P3/2 excited levels are 
separated by an angular frequency u{. The S1/2 electronic 
ground state consists of two hyperfine levels F — I — 1/2 and 
F = I + 1/2. The relative energies of these two levels de- 
pends on the sign of the hyperfine constant Am and can vary 
between ion species (in this figure, Am is negative). The qubit 
is encoded in the pair of = states of the two F manifolds 
separated by an angular frequency ujq. Coherent manipula- 
tions of the qubit levels are performed with a pair of laser 
beams that are detuned by A from the transition to the P1/2 
level, represented by the two straight arrows. The angular fre- 
quency difference between the two beams equals the angular 
frequency separation between the qubit levels uob — <^r = 1^0 ■ 
Some ion species have D levels with energies below the P man- 
ifold. Wavy arrows illustrate examples of Raman scattering 
events. 



Gates are assumed to be driven by pairs of Raman 
beams detuned by A from the transition between the 
S1/2 and the P1/2 levels (See Fig. [T]). We further assume 
that the Raman beams are linearly polarized and Raman 
transitions are driven by both a + photon pairs and cr_ 
photon pairs. The two beams in a Raman pair are des- 
ignated as red Raman (r) and blue Raman (b) by their 
respective frequencies. In the following we also assume 
that A is much larger than the hyperfine and Zccman 
splitting between levels in the ground and excited states. 

The Rabi frequency between the two clock states is 23 1 
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where g b/r = E b/r \(P 3/2 ,F = I + 3/2,m F = F|d • 
<7 + \St/2,F = I + 1/2, mj? = F)\/2H, E b/r is the peak 
electric field of the b or r beam at the position of the ion, 
respectively, and d • <7 + is the right circular component of 



the electric dipole operator. The right and left circular 
polarization components of the b and r beams are 
and r + /_, respectively. The Pi/ 2 and P 3 / 2 excited levels 
are separated by an angular frequency cjf. 

The total spontaneous photon scattering rate from 
these beams is given by [23| 



^ = l[9l(bl+bl) +9 l{rl+rl)\ 
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Here 7 is the natural linewidth of the P1/2 and P3/2 levels 
[I?]]. We now assume linearly polarized Raman beams 
with 6_ = b+ = r- = -r+ = 1/V2 We further 

assume g b = g r = g- The time for a 7r rotation is 
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Combining Eq. fl]),© and ©, the probability to scatter 
a photon during is given by 



"total — 



7T7 2A 2 + (A-u f ) 2 
W |A(A-Wf)| 



(4) 



The dashed line in Fig. [5] shows P to tai vs. A, where the 
laser detuning is expressed in units of the excited state 
fine-structure splitting, and the scattering probability is 
given in units of 7/wf. The total scattering probability 
has a global minimum of 



Pmin = 2V27T7/w f , 



(5) 



when the laser detuning is between the two fine-structure 
manifolds (A = (v2 — l)wf). The asymptotic value of 
^totai for large positive or negative detuning 



37T7/Wf, 



(6) 



is only slightly larger than the global minimum. 

Previous studies estimated decoherence by assuming 
that any photon scatter ing will immediately decohere a 
hyperfine superposition [23j . Under this assumption the 
lowest possible gate error equals P m in and ions with a 
small 7/wf ratio benefit from a lower gate error minimum. 

Two kinds of off-resonance photon scattering occur in 
the presence of multiple ground states. Inelastic Raman 
scattering which transfers population between ground 
states, and Rayleigh elastic scattering which does not 
change ground state populations. Since no energy or 
angular momentum are exchanged between the photon 
and the ion internal degrees of freedom, no information 
about the qubit state is carried away by a Rayleigh scat- 
tered photon. Rayleigh scattering, therefore, does not 
necessarily lead to decoherence [2J, [U [25[. This was 
experimentally shown in [25| where, when of equal rate 
from both qubit levels, off-resonance Rayleigh scattering 
of photons did not affect the coherence of a hyperfine 
superposition. Decoherence in the presence of light was 
shown to be dominated by Raman scattering. The guide- 
line used in [23j is therefore overly pessimistic. 
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FIG. 2: The solid line is the probability (Pse = ftuman)to 
scatter a Raman photon (in units of y/u)f) during a n ro- 
tation vs. the laser detuning in units of the excited state 
fine-structure splitting. The dashed line is the probability 
for any type of scattering event (Pse = Ptotai) during the 
pulse vs. detuning. The Raman scattering probability decays 
quadratically with A for |A| ^> cot. 



In this paper we re-examine the errors due to sponta- 
neous photon scattering on single-qubit gates (rotations 
of the equivalent spin 1/2 vector on the Bloch sphere) 
and two-qubit (entanglement) gates, where the qubits 
are based on ground state hyperfine levels and manipu- 
lated with stimulated Raman transitions. We compare 
between different ion species and examine different Ra- 
man laser parameters. In section [II] we analyze the con- 
tribution to the gate error due to spontaneous Raman 
scattering. Following a Raman scattering event the ion- 
qubit is projected into one of its ground states and spin 
coherence is lost. This error was also addressed in [26| . 
In section Hill we examine the error due to Rayleigh scat- 
tering. Rayleigh scattering error results primarily from 
the photon's recoil momentum kick. We show that both 
types of gate errors can be reduced to small values, while 
keeping the gate speed constant, with the use of higher 
laser intensity. 



II. RAMAN SCATTERING ERROR 

In a Raman photon scattering event energy and an- 
gular momentum are exchanged between the scattered 
photon and the ion's internal degrees of freedom. The 
polarization and frequency of the scattered photon (with 
respect to those of the laser) become entangled with the 
ions' internal state. Therefore, after tracing out the pho- 
ton degrees of freedom, the ions' spin coherence is lost. In 
other words, Raman scattering serves as a measurement 



of the ion-qubit internal state. Following a spontaneous 
Raman scattering event the ion-qubit is projected into 
one of the ground states in the S1/2 manifold. For ions 
with low lying D levels, Raman scattering events can also 
transfer the ion from the qubit levels into one of the D 
levels. 



A. Single-qubit gate 

For a single qubit gate we choose to look at the fidelity 
of a 7T rotation around the x-axis of the Bloch sphere, 
represented by the Pauli operator a x . This gate is as- 
sumed to be driven by a co-propagating Raman beam 
pair, where the frequency difference between the beams 
equals the frequency separation between the two qubit 
states loq. It is straightforward to generalize this case to 
other rotations. 

As a measure of the error in the rotation, we use the 
fidelity of the final state (characterized by density matrix 
Pfinai) produced by the erroneous gate as defined by 



F= (*|ftaal|*>, 



(7) 



where \*S>) is the ideal final state. Given an initial state 
l^init), I*) can be written as 



I*) 



(8) 



In the presence of off-resonance Raman scattering, the 
density matrix that describes the state of the qubit after 
the gate has the form 



Pfinal = (1 - fkaman)|*)(#| + Pe 



(9) 



The erroneous part of the density matrix p t — 
J2i i s composed of projectors into different levels 

\i) of lower energy. Here PR aman = J2i w i is tne proba- 
bility for a spontaneous Raman scattering event to occur 
during the gate. 

Note that some Raman scattering events keep the ion 
within the qubit manifold. Using Eq. ([7]) the contribu- 
tion of p e to the fidelity is positive and not strictly zero. 
For simplicity, we neglect this contribution and put a 
lower bound on the gate fidelity 



F ^ 1 pRaman • 



(10) 



In what follows we assume this expression to be an equal- 
ity. The error in the gate due to spontaneous Raman 
photon scattering is hence given by the Raman scatter- 
ing probability 



1 - F = P, 



Raman ■ 



(ii) 



We first examine the error due to Raman scattering 
back into the 2 Si/ 2 manifold eg. The Raman scattering 
rate back into the S1/2 manifold is calculated to be [25[ 



TRaman = y + b 2 + ) + + r\ )] 
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For the same laser parameters as above, the probability 
to scatter a Raman photon during the gate is 



Pi 



Raman 



2-7T7 



3 |A(A-w f )| 



(13) 



The solid line in Fig. [2] shows Pnaman vs. A. The 
Raman scattering probability decays quadratically with 
A for |A| 3> u>f. Qualitatively, this is because Raman 
scattering involves a rotation of the electron spin. Elec- 
tron spin rotations are achieved through the spin-orbit 
coupling in the excited state. This coupling has oppo- 
site sign contributions from the two fine-structure levels. 
Therefore as we detune far compared to the fine structure 
splitting, those two contributions nearly cancel (24[. 
Using Eq. (JTJ) we can write 



Pi 



Raman 



(14) 
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Gate Error 



The ratio can be expressed in terms of atomic con- 
stants, and the peak electric field amplitude E at the 
position of the ion [27j 



itkJ 3 3/2 



g 2 3neoc 3 E 2 



(15) 



Here w 3 / 2 is the frequency of the transition between the 
Si/2 an d the P3/2 levels (D2 line), eo is the vacuum per- 
mittivity, and c is the speed of light. Assuming Gaussian 
laser beams, at the center of the beam 



E 2 = 



4P 



7rwgce 



(16) 



Here V is the power in each of the Raman beams and wo 
is the beam waist at the position of the ion. The prob- 
ability to scatter a single Raman photon can be written 
as 



Raman 



2Tr\S}-n\h(jJ 3 ,/ 2 WQ 
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(17) 



This result is essentially the same as Eq. (5) of pq . 
We can now rearrange this expression to put an upper 
bound on the required power for a desired gate speed 
and error 



V = — 



2ir 2-kwq , 



3es A3/2 



' 7^3/2 1 I 



(18) 



where A 3 / 2 = c/uj 3 / 2 . 

Assume that the ratio of the beam waist to the tran- 
sition wavelength is constant for different ion species. In 
this case, the power needed to obtain a given Rabi fre- 
quency and to keep the error below a given value would 
scale linearly with the optical transition frequency. A 
more realistic assumption might be that the Raman beam 
waist is not diffraction limited and is determined by other 
experimental considerations, such as the inter-ion dis- 
tance in the trap or beam pointing fluctuations. In this 



FIG. 3: Laser power in each of the Raman beams vs. the error 
in a single ion gate (n rotation) due to Raman scattering back 
into the S1/2 manifold (obtained using Eq. (JTSJ). Different 
lines correspond to different ion species (see legend). Here 
we assume Gaussian beams with wo = 20 fim, and a Rabi 
frequency Q,-[{j2n = 0.25 MHz (tv = 1 fisec). 



case, assuming that wq is constant, the required power 
would scale as the optical transition frequency cubed. 
Either way, ion species with optical transitions of longer 
wavelength are better suited in the sense that less power 
is required for the same gate speed and error require- 
ments. In addition, high laser power is typically more 
readily available at longer wavelengths. Finally, we note 
that the error is independent of the fine-structure split- 
ting as long as we have sufficient power to drive the tran- 
sition. The transition wavelengths of different ions are 
listed in Table |U 

Figure [3] shows the laser power needed per Raman 
beam for a given error due to spontaneous Raman scat- 
tering into the S1/2 manifold. Here we assume Q-r_/2tt = 
0.25 MHz (^ = 1 /usee), and wq = 20 /im. Different lines 
correspond to the different ion species listed on the figure 
legend. Table |TT] lists the error in a single ion-qubit gate 
due to Raman scattering back into the S1/2 manifold for 
the same parameters as in Fig. [3J and assuming 10 mW 
in each of the Raman beams. The power Vq needed in 
each of the gate beams for es = 10~ 4 is also listed in the 
table. As can be seen, ions with shorter transition wave- 
length require a larger detuning and accordingly higher 
laser power to maintain a low gate error. For most ions, 
a few milliwatts of laser power is enough to reduce the 
gate error to below 10~ 4 . 

Equation (flU)) can be solved to give the required detun- 
ing values for a given es ■ The number of such detuning 
values comes from the number of crossings of a horizontal 
line, set at the desired error level, with the solid curve 
in Fig. [21 When es is higher than the minimum Pa am an 
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TABLE I: A list of atomic constants of several of the ions considered for quantum information processing. Here I is the nuclear 
spin, 7 is the natural linewidth of the P1/2 level 28, 29, 30l |3l l [3z | . [33 L [34l l35l . ujq is the frequency separation between the two 
qubit states set by the hyperfme splitting of the S1/2 level [3a. [371 13a. l39l l40l. I41L l4zl |4^. El] . u>t is the fine-structure splitting 
HH, A1/2 an d A3/2 are the wavelengths of the transitions between the S1/2 and the P 1/ 2 and P3/2 levels [45l ]. respectively. The 
branching ratio of decay from the P levels to the D and the S levels is / [4r3 . [47l |48| . |49| . [50l . [5ll . |52| . 



inside the fine-structure manifold, i.e. es > 87T7/3u;f (see 
Fig. [J) , then four different detuning values yield the same 
es, two outside and two inside the fine-structure mani- 
fold. When eg is lower than this value, only two detuning 
values, both of which are outside the fine-structure man- 
ifold, yield £3. Those detuning values An that are below 
the S1/2 — > P1/2 transition (A < 0) and correspond to 
es = 10~ 4 are listed in Table |TT] for different ions. For 
most ions An is in the few hundred gigahertz range, and 
its magnitude is much smaller than uif. 

We now consider the error for various ions caused by 
Raman scattering into low lying D levels, en- As tran- 
sitions between levels in the S and the D manifolds do 
not necessarily involve electron spin rotations, the error 
suppression discussed preceding Eq. (flU]) will not occur. 
Instead, Raman scattering rate into the D levels will be 
given by the total scattering rate times a fixed branch- 
ing ratio /. When driven resonantly the P1/2 and P3/2 
levels decay to the D manifold with different (but often 
similar) branching ratios. Here we assume that for a de- 
tuning large compared to the fine-structure splitting, the 
branching ratio is essentially independent of the laser de- 
tuning and is given by the average of the two resonant 
branching ratios [4(1 S3- Table [J lists/, for various ion 
species, obtained from j48l.l49l.[50l[5lT[52|. The error due 
to Raman scattering into D levels is given by 

e D = /Ptotal- (19) 

Using Eq. fUS]) we can write the ratio of the errors due 
to Raman scattering into the different manifolds 

For |A| < u>{, Raman scattering error is dominated by 
scattering back into the Si/ 2 levels. For ions with 1// >> 
1 the two errors become comparable at a detuning A ~ 
V / zV | f/3v'X When the detuning becomes large compared 
to the fine-structure splitting, scattering into low-lying D 
levels dominates. For most ions considered here, |Aq| < 
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TABLE II: A list of errors in a single-qubit gate (ir rotation) 
due to spontaneous photon scattering. The error due to Ra- 
man scattering back into the S1/2 manifold es is calculated 
with the same parameters as Fig. [3] Gaussian beams with 
Wo = 20 /im, a single ion Rabi frequency Qr/2tt = 0.25 MHz 
[t-it — 1 /isec), and 10 mW in each of the Raman beams. Vo 
is the power (in milliwatts) needed in each of the beams, and 
Ao/27r is the detuning (in gigahertz) for es = 10 -4 . The ra- 
tio between errors due to Raman scattering to the D and S 
manifolds, en/es, is given when es = 10 -4 . The asymptotic 
value of en in the jA| ^> u>{ limit is eooo- 



Wf and (perhaps with the exception of 137 Ba + ) eg is the 
more dominant source of error. The ratios eo/eg when 
es = 10~ 4 (i.e. A = An) are given in Table ITT1 for different 
ions. 

Due to the asymptotic value of the total scattering rate 
in the |A| ^> ojf limit (Eq. ©), cd has an asymptotic 
value which gives a lower bound to the Raman scattering 
error 

e Doo = ^. (21) 

Table [III lists cdoo for various ions. For all ion species 
considered this value is below the assumed estimates for 
the fault tolerance threshold. 



B. Two-qubit gate 

A universal quantum gate set is complete with the ad- 
dition of two-qubit entangling gates. During the last few 
years there have been several proposals and realizations 
of two ion-qubit gates Q, H, ©, ES EH El H EH- 
Here we focus on gates that use spin-dependent forces in 
order to imprint a geometric phase on certain collective 
spin states [1 [1 ElL El IH H3, HH • We examine only 
gates that are implemented with a continuous non reso- 
nant pulse rather than those using multiple short pulses 
[Ell ]. Again, to compare different ion species we examine 
ion-qubits that are encoded into hyperfine clock states. 
When the laser detuning is large compared to the hyper- 
fine splitting, the differential light force between clock 
levels is negligible @, However, a phase gate can be 
applied between spin states in the rotated basis (super- 
positions of clock states that lie on the equatorial plane 
of the Bloch sphere) [HI H3| ■ In this scheme the ions tra- 
verse a trajectory in phase space that is conditioned on 
their mutual spin state (in the rotated basis 11]). The 
phase the ions acquire is proportional to the total area 
encircled in phase space. This geometric phase gate was 
demonstrated in [8| and was realized on clock states in 

This form of phase gate is implemented with two dif- 
ferent Raman fields that are slightly off-resonance with 
upper and lower motional sidebands of the spin-flip tran- 
sition. For simplicity, we assume here that the gate is 
driven by two independent pairs of Raman beams, i.e. 
a total of four beams. Most experimental implementa- 
tions of this pha se gate thus far have used a three-beam 
geometry p| It is straightforward to generalize the 
treatment below to the three-beam case. 

Typical conditions for the gate are such that the angu- 
lar frequency difference between the beams is (wo +wt rap — 
6) in one Raman pair and (ujq — CL> trap + ^) m the other. 
Here Wtrap is the angular frequency of the normal mode 
that the gate excites and 6 is the Raman field detuning 
from that motional sideband. Under these conditions the 
ions will traverse K full circles in phase space for a gate 
duration of T ga t = 2ttK/S. Typically \S\ is chosen to be 
much smaller than uitmp to avoid coupling to the pure 
spin flip ("carrier") transition or the motional "specta- 
tor" mode, and K is usually chosen to be 1 to minimize 
the gate time. 

Figure [4] depicts the assumed geometry of the laser 
beams. The two beams comprising each Raman pair in- 
tersect at right angles at the position of the ions, such 
that the difference in their wave vectors is parallel to the 
trap axis. With this choice, the Raman fields couple, to a 
very high degree, only to the motion along the trap axis. 
The beam polarizations are assumed to be linear, perpen- 
dicular to each other and to the magnetic field axis. The 
beams' relative frequencies can be arranged such that 
the final state is insensitive to the optical phase at the 
ions' position [57| . Generalizing this treatment to other 
Raman beam geometries is straightforward. 




Laser beams 



FIG. 4: Schematic of Raman laser beam geometry assumed 
for the two-qubit phase gate. The gate is driven by two Ra- 
man fields, each generated by a Raman beam pair. Each pair 
consists of two perpendicular beams of different frequencies 
that intersect at the position of the ions such that the differ- 
ence in their wave vector lies parallel to the trap axis. One 
beam of each pair is parallel to the magnetic field which sets 
the quantization direction. The beams' polarizations in each 
pair are assumed to be linear, perpendicular to each other 
and to the magnetic field. Wavy arrows illustrate examples 
of photon scattering directions. 



As in the single-qubit case, Raman scattering will 
project one of the ion qubits into one of its states be- 
low the P manifold. In the appropriate basis the ideal 
gate operation is represented by [9|, lllj 
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where typically tj> = ir/2. The output state is |\&) = 
U\ ^init) • The density matrix following the erroneous gate 
will be of the form 



/Ofinal = (1 - 2P] 



Raman-gate 



|*><*| 



(23) 



where -PRaman-gato is the probability that one of the ions 
scattered a Raman photon during the gate [58|, l59( . For 
a gate consisting of K circles, the Raman detuning 6 is 
chosen such that the gate time is given by [oil ] 
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Here rj = Akzo is the Lamb-Dicke parameter, where 
Afc = is the wave vector difference between the 

two beams that drive the gate (for the particular geom- 
etry of Fig. 2]), and is the laser beam wave vector 



magnitude. The root mean-square of the spatial spread 
of the ground state wave function of one ion for the nor- 
mal mode that the gate excites is 



Z ° = \/ S-/4Mw t ra 



(25) 



where M is the mass of an individual ion. The single ion 
carrier Rabi frequency Or is given in Eq. |T]). 

Using similar considerations to those used in the single 
ion gate and, as before, neglecting any (positive) contri- 
bution of p e to the fidelity, we can set an upper bound 
for the power needed for a certain error in the gate due 
to Raman photon scattering into S1/2 levels 
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This required gate power is AyK/rj times larger than 
that needed for the same error in the single ion-qubit it 
rotation given in Eq. (JTSJ) . The factor of \fKjr\ is due to 
the longer two ion-qubit gate duration compared to single 
qubit rotations, and the factor of 4 is due to the presence 
of two ions and the pair of required Raman fields. 

When comparing different ion species we can fix differ- 
ent parameters, depending on experimental constraints 
or requirements. For example, here we choose as fixed 
parameters the beam waist wq (for the reasons given 
above), the gate time r ga t (assuming a certain compu- 
tation speed is desired), and the mode frequency Wtrap 
(which sets the time scale for various gates) . With these 
choices, heavier ions pay the price of smaller i] and there- 
fore higher power requirements per given gate time and 
error. A different approach would be to choose rj fixed, 
in which case heavier ions will need a lower Wtrap- For 
a fixed gate time a lower Wtrap leads to stronger, off- 
resonant, coupling of the Raman fields to the carrier or 
the other motional "spectator" mode, and, therefore, to 
a larger error due to this coupling [60, [fill . The Lamb- 
Dicke parameters for the different ions for wtrap/27r = 5 
MHz are listed in Table HTT1 

With the choice where wq, T ga t e , and Wtrap are fixed, 
we can write 
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Figure [5] shows the power needed vs. error in a two-qubit 
gate due to Raman scattering back into the S1/2 manifold 
for wo = 20 /im, T ga t c = 10 (is, wtrap/27r = 5 MHz and 
K = 1. Table IIIT1 lists es for the various ion species and 
the same laser parameters as Fig. [5j assuming a power 
of 10 mW is used in each of the four Raman beams. 
Alternatively, the power and the detuning A needed 
in each of the gate beams for es = 10~ 4 are also listed 
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An /PirfTH^ 
ljO / Z/l 1 1 11Z 1 


£d/ es 


fcDoo / -LU 


V 


a Be + 


3.6 


36 


-1.20 


N.A. 


N.A. 


0.194 


^ b Mg + 


11.1 


111 


-7.28 


N.A. 


N.A. 


0.130 


Ca + 


13.6 


136 


-10.42 


1.01 


1.06 


0.071 


Zn + 


41.1 


411 


-24.96 


N.A. 


N.A. 


0.11 


Sr + 


26.5 


265 


-20.34 


0.52 


0.50 


0.048 


m Cd+ 


64.3 


643 


-35.44 


N.A. 


N.A. 


0.081 


iyv Ba+ 


37.4 


374 


-30.67 


1.65 


1.46 


0.034 


m Yb+ 


57.5 


575 


-32.89 


0.01 


0.007 


0.038 


ly9 Hg+ 


149.7 


1497 


-49.52 


0.003 


0.001 


0.078 



TABLE III: A list of different errors in a two-qubit phase 
gate due to spontaneous photon scattering. The error due to 
Raman scattering back into the 1S1/2 manifold, es, is calcu- 
lated assuming Gaussian beams with wo=20 /im, a gate time 
Tg a tc=10 fis, £j t rap/27T =5 MHz, a single circle in phase space 
(K=l), and 10 mW in each of the four Raman beams. Vo 
is the power in milliwatts needed in each of the beams, and 
Ao/27r is the detuning in gigahertz for es = 10 -4 . The ratio 
between errors due to Raman scattering to the D and S man- 
ifolds, eo/es, is given when es = 10~ 4 . The asymptotic value 
of eo in the |A| S> uif limit is eDoo- The Lamb-Dicke param- 
eter rj for the above trap frequency is also listed for different 
ions. 



in the table. Here, heavier ions need a larger detuning 
and correspondingly higher laser power to maintain a low 
gate error. For most ion species, hundreds of milliwatts 
of laser power per beam and a detuning comparable or 
even larger than u>{ are needed to reduce the gate error 
to the 10~ 4 level. 

As in the one-ion gate, Raman scattering into low lying 
D levels will add to the gate error. This error is given by 



e D = 2/P t 



total-gate 



fPi 



total ■ 



(29) 



Here Ptotai-gatc is the probability that one of the ions 
scattered a photon during the two-qubit gate, and Ptotai 
is the one-qubit gate scattering probability given in Eq. 
(TJJ . Since both the Raman and the total scattering prob- 
abilities increase by the same factor as compared to the 
one-qubit gate, the ratio of the two errors eu/^s will re- 
main the same as given by Eq. (f20|) . Table HTI1 lists en/es 
for the different ions when es = 10 -4 . Notice that for es 
= 10 some ions require |Ao | > V2uj{/3^. For those 
ions cd is no longer negligible compared to es- 

Scattering into a low lying D level will, again, set a 
lower bound on the total error. Combining Equations 
([29)) and ([6]) we find this lower bound to be 



3tt 7 / WK 



ujf rj 

Table Hill lists eDoo for the different ion species. 
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FIG. 5: Laser power in each of the Raman beams vs. the error 
in a two ion entangling gate due to Raman scattering back 
into the S\/2 manifold (obtained using Eq. (|27|l 'l. Different 
lines correspond to different ion species (see legend). Here we 
assume Gaussian beams with u>o = 20 /im, a gate time r ga t e 
= 10 fis, aj t rap/27T = 5 MHz, and a single circle in phase space 
(K = 1). 

III. RAYLEIGH SCATTERING ERROR 

Since Rayleigh photon scattering is elastic, no en- 
ergy or angular momentum is transferred between the 
photons' and the ions' internal degrees of freedom. 
Therefore, these degrees of freedom remain uncorrelated. 
Rayleigh scattering does not necessarily lead to direct 
spin decoherence. 

In situations where Rayleigh scattering rates from the 
two ion-qubit states are different, Rayleigh scattering of 
photons will eventually measure the qubit state and lead 
to decoherence. In fact, the most common ion-qubit de- 
tection method relies on state selective Rayleigh scat- 
tering of photons on a cycling transition. In most Ra- 
man gates, however, the laser is detuned from resonance 
by much more than the qubit levels' energy separation, 
A 3> Wo, and the laser polarization is typically linear to 
suppress differential Stark shifts. Under these conditions 
Rayleigh scattering rates from the two qubit levels are 
almost identical and the error due to the rate difference 
is negligible. For the clock transition qubit states consid- 
ered here, Rayleigh scattering rates are almost identical, 
regardless of the laser polarization (for more details on 
this error see subsection IIIIBJ) . 

The main effect of Rayleigh scattering is the momen- 
tum recoil it imparts to the ion-qubit. For a single- 
qubit gate, the Raman beams are usually arranged 
in a co-propagating geometry. In this configuration, 
since the Lamb-Dicke parameter is very small (77 ~ 
(wo/w3/2)fcL<Zo), the effect of ion motion is negligible on 



the gate operation. Therefore, Rayleigh scattering has a 
negligible effect on single-qubit gates. 

In the two-qubit phase gate, a mode of motion is ex- 
cited that is entangled with the two-ion collective spin 
state. In this case, recoil from photon scattering perturbs 
the ion's motion through phase space and contributes to 
the gate error. The ion-qubit trajectory is distorted in 
two ways. The larger distortion arises from the direct re- 
coil momentum displacement. A second, much smaller, 
distortion arises from the contribution of the recoil to the 
appearance of nonlinearities in the gate operation due to 
deviations from the Lamb-Dicke regime (see subsection 

MM- 

In subsection IIII Al we calculate the error in a two- 
qubit gate due to direct recoil phase-space displacement. 
In subsection lHIBI we elaborate on the two other sources 
of error mentioned above, namely errors due to uneven 
Rayleigh scattering rates and errors due to deviations 
from the Lamb-Dicke regime. 



A. Rayleigh scattering recoil error 

In elastic Rayleigh scattering, energy and momentum 
are not exchanged between the ions' and the photons' 
internal degrees of freedom. However, momentum and 
energy are exchanged between the photons' and the ions' 
external degrees of freedom. The scattered photon direc- 
tion will be different from that of the laser beam, caus- 
ing the ion to recoil. This recoil acts on the ion-qubit 
as a phase-space momentum displacement, distorting the 
ion's trajectory through phase space and causing an error 
in the entangling-gate phase. Note that the momentum 
imparted to the ion (and therefore the deviation from 
the desired gate phase) and the momentum that is car- 
ried by the scattered photon (namely its scattering di- 
rection) are correlated. From this point of view the gate 
infidelity again arises due to the entanglement between 
the scattered photons' and ions' (this time external) de- 
grees of freedom (for more details on this point of view 
see Appendix [Af . 

In the Lamb-Dicke regime (for a thermal state 
i]\/2h + 1 <C 1, where n is the average mode popula- 
tion) , the gate operation can be approximated as a series 
of finite displacements in phase space 0] 

N 

U=l[D{Aa k ). (31) 
fc=i 

Displacements through phase space, Aa^, are condi- 
tioned on the joint spin state of the two ions, and depend 
on the gate parameters. For certain gate parameters, the 
displacement is zero for the two parallel spin states | tt) 
and I !!), and non zero with opposite sign for the | t-L) 
and I !f) states, where | and j hereafter refer to the ro- 
tated basis rather than the clock levels. Using the com- 



mutation relations between phase space displacements 

b{a)b{J3) = D(a + /?)e lIm(Q/r) , (32) 

we can write the gate operation as a single displacement 
times an overall phase 



U 



N 
k=l 



(33) 



where N is the total number of infinitesimal dis- 
placements. When the net displacement is zero, i.e. 
X^jtLi Aaj. = 0, the two-ion motion returns to its ini- 
tial state and spin and motion are disentangled at the 
end of the gate. The gate operation on the affected spin 
states can be written as 



The phase acquired, 



U = Ie 



N 



%4> 



= Im(^Aa J ^Aa / ), 



(34) 



(35) 



J=2 



7 = 1 



is the same for the | Tl) and | J.T) states and proportional 
to the encircled phase-space area. Figure^ shows an ex- 
ample for the trajectory traversed in phase space during 
an ideal gate in a reference frame rotating at the mode 
frequency w t ra P - 

In the presence of laser light there is a finite probability 
fkayieigh-gate that a Rayleigh photon will be scattered by 
one of the ions during the gate [58[ . 

The probability for a Rayleigh scattering event to oc- 
cur during a single-qubit gate is given by the difference 
between Eq. (0} and Eq. ^ 



Rayleigh 



3A 2 - 2Au f + ujf/3 



|A(A-o*)| 



(36) 



Since in the limit of |A| ^> u>{ all scattering events are 
Rayleigh scattering, PRayieigh and Ptotai have the same 
asymptotic value P^. 

Using Eq. (f36|) and the factor for the extra required 
power in the two-qubit gate 



Raylcigh-gatc 



AVKttj /3A 2 - 2Aw f + w f 2 /3 
Van \ |A(A- w f )| 



(37) 



The effect on the ion motion would be that of momentum 
recoil 



(38) 



Here q = Icl — kg cat is the wave vector difference between 
the scattered photon and the laser beam from which it 
was scattered, and fj is the position operator of the ion 
that scattered the photon (i = 1,2). We will neglect 
recoil into directions other than along the trap axis. We 




FIG. 6: Schematics of the trajectories traversed by the ions 
in phase-space during the gate. The acquired phase is pro- 
portional to the encircled area, (a) Phase-space trajectory 
of an ideal gate, (b) Phase-space trajectory of an erroneous 
gate where a photon was scattered during the gate drive. The 
short straight arrows represent the recoil displacement. The 
grey-shaded area is proportional to the phase error. The area 
that is added to the | jj.) trajectory is subtracted from the 
| |t) trajectory. 



can write the position operator for for ion i along this 
axis as 



Zi + z (a + aj) + zi(ai + a|), 



(39) 



where Zi is the equilibrium position of the ion and Zo/ii 
aj^ and a /i arc the root mean-square of the ground 
state spatial spread, and the creation and annihilation 
operators, respectively, of the two motional modes and 
1. We assume the gate is performed by exciting mode 0. 
Recoil into mode 1 does not distort the gate dynamics 
directly. This part of the photon recoil will add minutely 
to the gate error through its contribution to the gate 
nonlinearity discussed in section IIII Bl We neglect this 
contribution for the moment and write the recoil opera- 
tion as 



Un 



O iq z 2o(a +aj) _ /3a -/3*aJ 



DOS), (40) 



where q z is the projection of q along the trap axis, and 
f3 = iq z z . 

The recoil action on the two-ion crystal is, therefore, 
a momentum displacement (3 in phase space. The mag- 
nitude of (3 depends on the wave-vector difference be- 
tween the scattered photon and the beam from which 
it was scattered, and the projection of this momentum 
difference along the trap axis. Since the same recoil dis- 
placement is applied to all spin states, motion will still 
be disentangled from spin at the end of the gate. Er- 
rors are therefore due to the change in area encircled in 
phase space. Figure [6Jd illustrates a gate trajectory that 
is distorted due to photon recoil. 

The erroneous gate can be again written as a sum of 
displacements, those due to the gate drive and (3, which 
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occurs at some random time during the gate. For a par- 
ticular two-ion spin state, 



N 



M 



Ue= H D(Aa k )D(p) Y[D(A ai ). (41) 



k=M+l 



1=1 



By use of the commutation relation (|32|) , the gate is writ- 
ten as 

U t = D(f3y^ +A *\ (42) 
where the phase error A<p is determined to be 



A(j) = Im 



M 



N 



(43) 



1=1 k=M+l 

Writing displacements as a function of the gate time 



Aa k = Aa(t k ) -> ^pdt, 



(44) 



we can write the phase error as a function of the gate 
displacement before and after the scattering time £ scat 



A(j) = Im 



(45) 



dt ' J, dt 

Since a(t gatc ) = ot(0) — we get 

A(b = Im[/?a*(t scat ) - 0*a(t scat )]. (46) 



Since the accumulated displacement at the moment the 
photon was scattered a(i scat ) is of equal magnitude but 
opposite sign for the | f!) and | if) states, we get 



A(j) U = -A0 iT = Acj>. 



(47) 



The erroneous gate can be therefore represented by the 
operator 



/l 0^ 

e 4 ^ +A ^ 

e *(<i>-M) o 

\0 1/ 



(48) 



The final state fidelity following the erroneous gate de- 
pends on l^init)- For a general initial state 



l*init. 



Ml TT> + «| T4)+7l ll) + s\ 11), 



(49) 



and using Eq. ([7]), we can write for the gate fidelity 



- ||„|2 j_ e iA*i ,2 + e -iA* U |2 



W + W , (50) 



The most relevant fidelity for fault tolerance considera- 
tions is that averaged over all possible initial states. The 
fidelity we calculate here is that due to the worst-case 
input state, that is the input state that minimizes Eq. 
(f50|) . The worst-case fidelity is clearly smaller than the 
average fidelity and therefore gives a conservative esti- 
mate to the gate error [62]. The worst-case fidelity for 
this kind of error was calculated in [63|. For A<j>/<j) < 1 
it is 



F = cos 2 A<t>. 



(51) 



This minimal fidelity is the result of an input state with 
|k 2 | = |7 2 |=1/2. 

Photon scattering occurs in only a small fraction 
-Prtayieigh-gate of the gates. Averaging over all gates per- 
formed, we get 

Fgate = 1 - -PRaylcigh-gatc(l " <COS 2 A(j>)) = 1 - CR, (52) 

where the brackets correspond to an average over all pos- 
sible A(f>'s, i.e. due to different values of f3 and different 
scattering times. 

To perform the above average we need to explicitly 
write the different displacements. In the rotating frame 
the cumulati ve g ate displacement after a time t can be 
written as 0, |53| 



«(*) 



-iSt 



2y/K 



l]e 



i<i>i 



(53) 



where is the gate phase, determined by the optical 
phase difference between the gate beams at the ion's po- 
sition. In the rotating frame, the recoil displacement can 
be written as 



= |/3|e l " trapt . 



(54) 



Substituting equations (j54"|) and ([53")) into Eq. (|4"6")l , we 
can write the averaged term in the fidelity 



(cos 2 A(j>) 



2n 



|/3|=0 Jt=0 



S(\0\) cos 2 



[C0s((w t rap + S)t - $i) - COs(w t rapi ~ 



dtd\0\. 



(55) 



Here we assume that the probability of scattering at dif- ferent time intervals during the gate is uniform. The 
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probability distribution for different recoil displacement 
magnitude |/3| is given by which is determined 

by the geometry of the Raman beams with respect to 
the trap axis and the probability distribution of photon 
scattering directions. 

With the laser beam configuration assumed in Fig. ([4]), 
recoil due to photon absorption can be imparted only at 
45° to the trap axis, whereas recoil due to photon emis- 
sion can be imparted in any direction. The maximum al- 
lowed displacement following an absorption-emission cy- 
cle is therefore 



= k L z a (l + V2)/V2 



2 (1 

2 V 



(56) 



For proper gate operation we require 77 <C 1; thus we ex- 
pand Eq. (155j) in powers of |/3|, including terms to order 
|/3| 2 . Further, since typically 5/wtrap <C 1, we neglect 
terms proportional to 5/wtrap- With these approxima- 
tions the gate fidelity is independent of $l and is given 
by 



E, 



ClgTl- 



gatc " 



IK 



(57) 



Here (|/3| 2 ) is the average of the recoil displacement mag- 
nitude squared 



2 ) 



|/3|: 



(58) 



For the assumed laser beam configuration, the light 
field polarization will oscillate rapidly during the gate 
between left /right circular and linear at the ions' posi- 
tion. Therefore, the scattering probability will average 
to be isotropic. For the beam configuration illustrated in 
FigHl in spherical-polar coordinates, 
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87T J Q 
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71 



cos ( 
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eded<t>. (59) 



Here the z (6 = 0) direction is chosen parallel to the trap 
axis. We find 



Combining Eqs. (JS2J), fl5 
error can be written as 



(60) 

and ([57]) . the Raylcigh 
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(61) 



(62) 



Table IIVI lists er/es f° r different ion species for 
w trap /27r=5 MHz, and a single-circle (K=l) gate, at a 
detuning where es — 10~ 4 . For most ions considered 



Ion 


fvi /fa 


ETi^/10 -4 


r>e 


1 A AO 

1. 442 


1 C 1 


Mg 


n 1 a o 
U.142 


U.104 


Oa 


0.0168 


0.018/ 




n m to 
0.01 1 2 


n m no 
O.Olyo 


8Vn + 

or 


U.UUoU 


U.UUo4 


m Cd+ 


0.0040 


0.0044 


iMV Ba+ 


0.0010 


0.0011 


ivi Yb + 


0.0006 


0.0006 


19y Hg+ 


0.0015 


0.0012 



TABLE IV: The ratio er/es of errors due to Rayleigh scat- 
tering recoil and Raman scattering to the S manifold, is given 
when es = 10 -4 . The asymptotic value of en in the |A| 3> Wf 
limit is ep,oo. Both are calculated assuming ai trap /27r = 5 MHz 
and K = 1. 



here, the gate error is still dominated by Raman scatter- 
ing at this level. 

For laser beam detuning large compared to the excited 
state fine-structure splitting (A 3> wt) scattering is dom- 
inated by Rayleigh events. In this limit, the probability 
that a Rayleigh photon will be scattered during the gate 
is 



Pi 



Raylcigh-gateoo 



v 



(63) 



The recoil error in a two-qubit gate is therefore asymp- 
totically bound by 



£Roo 



57T77 7 57T 2 7 



2VK v{ ojfX 



3/2 



2Mu tIav K 



(64) 



Table IIVI lists eRoo for different ion species for 
<Arap/27r=5 MHz and a single-circle (K=l) gate. With 
the exception of 9 Be + , the error due to photon recoil in 
a two ion-qubit gate is below 10 -4 . For this error heavier 
ions benefit due to their smaller recoil. 

It is possible to reduce £r by choosing a smaller 
rj/^J~K (by increasing the trap frequency and/or perform- 
ing multiple-circles gates). Correspondingly more laser 
power, proportional to K/i] 2 , is then required in order 
not to increase the Raman scattering error es or reduce 
the gate speed. 



B. Other Rayleigh scattering errors 



In section lTlI Al we calculated the error due to the recoil 
imparted to the ion-qubit during Rayleigh photon scat- 
tering. As noted, Rayleigh scattering of photons adds 
two more contributions to the gate error. 

The first contribution eg is due to a difference in the 
Rayleigh scattering rates from the two ground state lev- 
els. Assume that the Rayleigh scattering rate is rR ay i e ; g h 
from one qubit level and rR ay i e i g h + dT from the other. 
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A measurement of the qubit level will be conclusive once 
the difference in the number of photons that are scat- 
tered is larger than the standard deviation of the number 
of photons scattered from each level. Since the number 
of photons that are scattered follows a Poisson distribu- 
tion, a measurement will occur after photons have been 
scattered for a time t such that 

dTt > ^/T Ray i cigh t. (65) 

The error rate is, therefore, given by dr 2 /TR ay ieigh, and 
the error during a gate by 

£<5 = (dr/rR ay l e i g h) 2 PRayloigh- (66) 

For the above gate parameters, the difference in scat- 
tering rates is due solely to the difference in detuning 
between the two qubit levels. For A ^S> u>o, which 
is typically required to reduce other scattering errors 
to low levels, this error can be approximated by eg ~ 
(wo/A) 2 i-Rayieigh- The contribution of e$ to the total 
error at a realistic laser detuning is very small. We per- 
formed an accurate calculation of the difference in the 
scattering rates between the two clock levels, and veri- 
fied that indeed, for all ion species and laser detunings 
discussed above, es is negligible. 

A second source of Rayleigh scattering error is through 
the contribution of the recoil momentum displacement to 
non-linearities in the gate evolution. The gate error due 
to non-linearities was calculated in [o"H and in the present 
context is proportional to i"Rayieigh-gatef? 4 Var(n), where 
Var(n) is the variance of the motional modes distribution 
due to recoil. The recoil momentum displacement mag- 
nitude is of order r\. Therefore, starting from the ground 
state and following a single scattering event, Var(n) ~ r\. 
The gate error due to this effect will be proportional to 
-Prtayieigh-gateT? 5 oc 7/ 4 , and significantly smaller than other 
scattering errors discussed above. 

IV. DISCUSSION 

We have calculated the errors due to photon scattering 
in a single-pulse single-qubit gates and two-qubit phase- 
gates implemented with stimulated Raman transitions. 
These errors present a fundamental limit to the gate 
fidelity in trapped-ion QIP experiments that use these 
kinds of gates, and should be a significant factor when 
choosing a specific ion as a quantum information carrier 
for fault-tolerant quantum computing schemes. 

Three main errors occur from spontaneously scatter- 
ing photons during a gate. Two errors are due to Ra- 
man scattering either back into the S1/2 manifold (es) 
or to low- lying D levels (eu)- The third error is due to 
the Rayleigh scattering recoil during a two-qubit gate 
(cr). For most ions currently considered for QIP experi- 
ments, the dominant error for realistic laser parameters is 
es- This error can be typically reduced to below current 
estimates for the fault-tolerance threshold with the use 



of relatively high (but probably attainable) laser power. 
This makes the availability of high power laser sources 
at the relevant wavelength important. This error is also 
reduced for ions with a relatively large S — > P transition 
wavelength. 

Among those three errors, only tu cannot be reduced 
below a certain value by the use of higher laser intensity. 
Therefore, it may eventually be advantageous to choose 
ions that do not have low lying D levels. However, for 
most ions considered here, £d is still small at realistic 
laser parameters compared to the other Raman error es- 

The masses of the different ions play an interesting role 
in the scattering error. Since it is harder to transfer mo- 
mentum to heavier ions, they suffer from a larger eg in 
two-qubit gates (or alternatively from the need for more 
laser power for a given gate speed and error level). For 
the same reason, lighter ions suffer from a larger Rayleigh 
recoil error cr. Examining Eqs. (|28p and (|64p. the Ra- 
man error scales as tUg-M^trap and the Rayleigh error 
scales as l/wfy^Mo/trap (neglecting differences in wave- 
length and natural linewidths). Since both wq and u>{ 
are generally larger for heavier ions, the Raman error is 
larger and Rayleigh error is smaller for heavier ions by 
more than is indicated by the kinetic argument above. 
Since for most ions es is the dominant error, lighter ions 
seem to currently have a lower overall error due to pho- 
ton scattering. Also, because both errors are a func- 
tion of -M^trap lighter ions will reach the same error level 
and power requirements with a higher entrap when com- 
pared to heavier ions. This allows for faster gate opera- 
tion. For a given trap geometry and applied potentials, 
the axial trap frequency scales as 1/yM and the radial 
trap frequencies scale as 1/M. For the ion crystal to re- 
main along the axial trap direction the radial frequencies 
have to be larger than the axial frequencies; therefore the 
limiting frequency is the radial frequency. In this case, 
Mwtrap is independent of the ions' mass. Heavier ions 
under these conditions will have a lower wtrap, leading to 
slower gate operation. 

In conclusion, to minimize the effect of scattering on 
the fidelity of trapped- ion-qubit gates, one needs to strike 
a balance between the desirable characteristics of long 
wavelength, light mass, the availability of high power 
laser sources and, if possible, the lack of low lying D levels 
when choosing a specific ion as a quantum information 
carrier. 

Finally, we remind the reader that here we have fo- 
cused on hyperfine ion-qubits and gates that rely on 
off-resonant Raman transitions applied in a continuous 
pulse. Other kinds of trapped-ion gates or qubits could 
have different limitations on the gate fidelity due to spon- 
taneous photon scattering. 

We thank E. Knill for many insightful discussions and 
comments, J. Home, J. Amini and A. S0rensen for helpful 
comments. This work is supported by DTO and NIST. 
Contribution of NIST; not subject to US copyright. 
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APPENDIX A: DEPHASING VS. 
ENTANGLEMENT WITH DIFFERENT PHOTON 
MODES 

In section [TTT] we calculated the error due to Rayleigh 
photon scattering without considering the scattered pho- 
ton degrees of freedom. Rather, we considered dephasing 
due to the random phase that is generated by the scat- 
tered photon's recoil. 

As noted, we can also view the gate error caused by 
Rayleigh scattering of photons as arising from the entan- 
glement between the photon and the ion-qubit external 
degrees of freedom. For example, a photon that is scat- 
tered during the gate, while the two spin wave packets 
are displaced from each other, can be collected by an 
ideal imaging system. Since the two-ion collective spin is 
entangled with their position, the position at which the 
photon is detected at the image plane can give "which- 
way" information, thereby measuring the ions' spin, col- 
lapsing the entanglement, and causing an error in the 
gate. Experiments investigating atomic decoherence due 
to this effect were performed in neutral-atom interferom- 
eters [HI, [65[ . This effect was also calculated in [66[ for 
a single ion in a superposition of two different motional 
coherent states. 

The general equivalence of these two points of view 
(dephasing vs. entangling with the environment) was ex- 
plained in [67} . Here we show this equivalence for the 
trapped ion case where the ion is undergoing an (ideally) 
spin-dependent closed loop displacement. For simplic- 
ity we examine the case of a single ion; it is, however, 
straightforward to generalize the following proof to the 
case of two ions. 

We first calculate the final state fidelity by examining 
entanglement with the photon modes. As in the two- 
qubit gate, a single ion is prepared in an equal superpo- 
sition of spin states and is cooled to the motional ground 
state: 

l*> = (^IT> + ^li»®|0>M||0) k . (Al) 

Here kets with subscript M represent motional states, 
and kets with the subscript k the different photon modes, 
which are initially empty (neglecting the laser mode). 
The ion is subsequently driven by an oscillating force 
which is detuned from its motional resonance. The direc- 
tion of force is opposite for the two different spin states. 
Ideally the two spin states would traverse opposite circu- 
lar trajectories in phase space (in a frame rotating at the 
motional mode frequency). At the end of the gate drive, 
both parts of the superposition (ideally) return to the 
ground state of motion, and acquire the same geometric 
phase 0gate which is proportional to the phase-space area 
encircled. The state at the end of the ideal gate drive is 

|*idcai) = e^(-^=| T> + -j=\ J.)) ® |0) M I |0) k . (A2) 

Now assume the ion scattered a photon during the gate 



drive. Immediately before the scattering event, the gate 
evolution produces 

l*> = (^| T>|«> + ±=e«\ |>| - a)) ® |0) k> (A3) 

where ±a are the conditional phase-space displacement 
for the two spin states (we have dropped the subscript 
M) and <fi is the geometric phase accumulated by the 
time of scattering. Immediately after scattering, a single 
photon is created in mode k' and the ion correspondingly 
recoils: 



I*) 



(A4) 

i- J 5 (^ k ,)(i= e ^e^|T)|a) + -^e^ , - f | |)|-a)) 
<8>|l)k> ® |0) k dk'. 

Here q' = k' — ki, is the wave vector difference between 
the scattered photon and the laser beam. The delta- 
function S(Ek') enforces energy conservation. The elec- 
tromagnetic field is now represented by a superposition of 
states in which a single photon was scattered into a cer- 
tain mode while all other modes are empty. Each part 
of this superposition is correlated with the corresponding 
momentum recoil operator acting on the trapped ion. Re- 
coil into directions other than the trap axis will give rise 
to motion which is common to both parts of the super- 
position, and can therefore be traced over. We, there- 
fore, neglect recoil into the dimensions other than the 



trap axis direction, and approximate e lq 



where 



q' z = q' • z. By use of Eq. (|40j) this momentum recoil can 
be written as a phase space displacement e 4q ' z = D(/3'), 
where (3' — iq' z zo, and zq is the root mean-square of the 
ground state spatial spread. Using Eq. ([3"!?]) we add the 
recoil to the gate displacement and the state after scat- 
tering can be written as 



I*) = 



(A5) 



i- J 8{E k >)(-j=e i +e <!3a W a ">\ T)|a + /3') 

+ -^=e^e- lIm ( /3 ' Q *)| |)| - a + /?')) ® |k')dk\ 

Since the total gate displacement is equal to zero, the 
remaining part of the gate displacement is =F« (depending 
of the spin state). Using Eq. (|3"!i|) again, we can write 
the state of the system after the gate has completed 



l*> = 
1 

4tt 



(A6) 



4tt J V2 
v2 

We can now define the gate phase-error Acj>' = Im(a/3'" 
a*/3'), and write the corresponding density operator 
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P=^ty J J 5(E k ,)6(E^)(e-^'e^"\ T)(T | + e iA *V iA *"| [){[ | + (A7) 

I 



Since we have no information about the mode into which 
the photon was scattered, we reduce the above den- 

I 



Note that the coherences of the density matrix in Eq. 
(|A8|) are given by an average of e lA * . A large spread in 
Acf>' will "wash-out" coherence and leave an incoherent 
statistical mixture. Tracing over the motional degrees of 
freedom we are now ready to evaluate the fidelity with 
respect to the ideal gate output state (Eq. (|A2|) ). 

F=\ + \ <cos(2A0')) = (cos 2 A0') , (A9) 

where averaging is performed with respect to all of the 
photon modes: 

(cos 2 A(//) = — / 6(E k >) cos 2 A^'dk'. (A10) 
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